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CL|' Abstract 

' The pre-exponential factor in the probability of decay of a metastable vacuum is 

calculated for a generic (2+1) dimensional model in the limit of small difference e of 
^ ' the energy density between the metastable and the stable vacua. It is shown that 

I this factor is proportional to e^^^^ and that the power does not depend on details of 

the underlying field theory. The calculation is done by using the effective Lagrangian 
method for the relevant soft (Goldstone) degrees of freedom in the problem. Unlike 
in the (1+1) dimensional case, where the decay rate is completely determined by the 
parameters of the effective Lagrangian and is thus insensitive to the specific details of 
the underlying (microscopic) theory, in the considered here (2+1) dimensional case the 
pre-exponential factor is found up to a constant, which does depend on specifics of the 
underlying short-distance dynamics, but does not depend on the energy asymmetry 
parameter e. Thus the functional dependence of the decay rate on e is universally 
determined in the considered limit of small e. 



The problem of decay of a metastable vacuum state has attracted interest since long ago 
both in statistical physics[Tl |21 El El and in the relativistic field theoryjHl El 13 El- In the 
latter setting the metastable (false) vacuum is a state of one or more fields corresponding to a 
local rather than the global minimum of the energy density. Such state is stable under small 
quantum fluctuations of the field(s), however it can decay into a lower-energy vacuum state 
through large fluctuations described by the quantum tunneling. The mechanism of such 
transition [HIIZ] is quite similar to that in the first order phase transitions and is described by 
nucleation and subsequent expansion of bubbles of the lower (true) vacuum amidst the bulk 
of the metastable phase. In the process of the expansion the excess e of the energy density in 
the metastable vacuum (latent heat) is being converted into the energy of the (expanding) 
surface of the bubble. Clearly, the classical expansion of the bubble is possible only if the 
gain in the volume energy e ■ [Volume) compensates for the positive energy of the surface 
of the bubble /i ■ [Surface Area) with /i standing for the surface tension. The minimal 
configuration where these are exactly equal is the critical spherical bubble with the radius 
Rc = [d — l)/i/e, where d is the total number of the space and time dimensions. Starting 
from Rc the bubbles expand classically. The evolution in the classically forbidden domain 
at R < Rc is described by the quantum tunneling. 

It should be noted that in the above description the thickness of the surface is totally 
ignored in comparison with the radius of the bubble. The scale for the thickness of the 
transition region between the phases is set by the Compton wavelengths of the relevant 
particles in either of the vacua. Assuming that no relevant massless particles in either of the 
vacua are present in this problem, one readily concludes that the thin wall approximation is 
always valid in the limit of small e which is assumed throughout the present paper. 

The rate F of nucleation of a critical bubble per unit time and per unit volume is deter- 
mined by the tunneling part of its trajectory in the classically forbidden domain, and the 
exponential factor, found from the action on this trajectory, is well known[Sl III in terms of 
fi and e: 

T = Fe-\ (1) 

where B = Vd [d—l)'^~^ /i'^/e'^"^ with Vd = ir'^^'^ /T[l + d/2) being the volume of a unit ball in 
d dimensions. A calculation of the pre-exponential factor F requires a summation over the 
fluctuations of the field(s) near the tunneling trajectory[H]. In a (1+1) dimensional model, 
i.e. for d = 2, such summation can be done entirely within the thin wall approximation, and 
the factor F in eq.(^ is given[5] hj F = e/(27r). The latter expression is universal in the 
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sense that the only quantities that determine the rate T are /i and e. In other words, any 
possible complexity of the underlying field dynamics reduces for the purpose of calculation 
of the rate F to two macroscopic parameters^ /i and e. The (1+1) dimensional expression 
for F fully agrees with calculations in specific field models pm ITTl IT^ . 

The purpose of the present paper is to consider the problem of calculating the pre- 
exponential factor F in a (2+1) dimensional theory, i.e. for d = 3. Unlike in the (1+1) 
dimensional case, the knowledge of only the macroscopic parameters and e is insufficient 
for a complete calculation of F for d = 3. However it will be shown that the dependence of 
the factor F on the parameter e in the limit e — > is still universal, while the dependence 
on (a combination of) the parameter fi and any mass scales in the underlying field theory, 
does depend on the details of the model (i.e. on the 'microscopic' dynamics). Thus some 
universality still remains in the (2+1) dimensional case, although in a substantially reduced 
form. Namely, it will be shown here that in (2+1) dimensions the rate of the false vacuum 
decay behaves at e — * as 



coupling constants, of the underlying field model, but does not depend on e. The parameter 
H in eq.© is the renormalized surface tension of the boundary separating the two vacuum 
phases in the limit e — 0. It can be noted that the power of e in the pre-exponent in eq.© 
is in agreement with the result of a direct calculation [13j in a specific 0^ model. The claim 
in the present paper is that this power is universal in the limit e so that it does not 
depend on details of specific model. 

The calculation of the tunneling exponential factor as well as of the pre-exponent is 
conveniently done in terms of the Euclidean-space formulation of the model[71 |H]. The 
partition function Z for the metastable vacuum state is written in terms of the Euclidean 
action S in the standard form of the path integral: 



where A/" is the normalization factor and the path integration runs over all the fields in the 

model, generically denoted here as 0, with the boundary condition that the fields approach 

^It can be also mentioned that the resulting equation in (1+1) dimensions has no corrections in powers 
of the dimensionless parameter e/^^. Only terms with higher exponents of — 7r/i^/e are possiblejHI- 
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where the dimensional constant A depends on the details, such as mass parameters and 
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their values in the metastable vacuum at the boundaries of the normahzation space-time 
box. The metastabihty of the considered 'vacuum' state results in that the partition function 
develops imaginary part (similar to the imaginary part of the energy of a resonance) which 
determines the decay rate (per unit space-time volume) according to the relation fOl [7] 
r = 2Im(ln Z) I (V^T) with VT being the space-time volume of the normalization box. 

The action functional S'[0] has a saddle point configuration corresponding to the tunneling 
trajectory of the bubble, the so-called bounce[Zj. In the thin wall approximation the bounce 
is a three-dimensional ball of the 'true' vacuum separated from the bulk of the false vacuum 
by the thin wall with the surface tension /i. The radius Rc of the bounce is thus determined 
by the extremum of the effective action 

S,jf = ^AB-tVB, (4) 

where Vb is the volume of the ball and Ab is its surface area. Clearly, the extremum of the 
action is achieved at Rc = 2yu/e, which is the critical bubble size in (2+1) dimensions, and 
the value of the action on the bounce is Sb = B = (167r/3) fi^/e^, which gives the exponential 
factor in eq.Q- 

In order to calculate the pre-exponential factor one has to evaluate the path integral Zi 
around the bounce configuration. The one- loop expression for the rate then reads as|HI 



VT 



det{S 



(2). -1/2 
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where Si is the operator of the second variation of the (Euclidean) action at the bounce 

(2) 

configuration, and 5*0 is the same operator for variation of the fields around the 'flat' false 

(2) 

vacuum state. The operator SI has three translational zero modes and the integration over 
those cancels against the space-time normalization factor in the denominator. In Ref . ^ this 
integration has been done explicitly. However in the present discussion it is somewhat more 
convenient to leave the expression in the symbolic form (jSJ and to deal with the zero modes 

(2) 

later. Also it should be mentioned that the operator Si has one negative mode, integration 
over which produces the imaginary part of the partition function, as explained in Ref.|H]. 

The spectrum of the operator S^ consists of two substantially different parts: the hard 
part with the eigenvalues A„ starting at the scale set by the mass parameters of the underlying 
field model, which scale is generically denoted here as m, and the soft part, whose scale is 
set by the (inverse) radius of the bounce. The soft part of the spectrum is universal and 
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does not depend on details of the underlying field model as long as the condition mRc ^ 1 is 
satisfied, which is always the case in the thin wall limit, i.e. at e — > 0. The explicit expression 
for the soft eigenvalues is[S] 

A.= , (6) 

where £ = 0,1, 2,. ..is the angular momentum and, obviously, the degeneracy of each mode 
is 2£ + 1. It can be readily noticed that the soft eigenmodes coincide (up to an overall nor- 
malization) with those found from the effective action Indeed, using the parametrization 
of the surface of the bounce in polar coordinates as r{6, (f) and expanding the radial position 
r around its value at the extremum of the action: r{6,ip) = Rc + a{6,ip), one gets in the 
quadratic order in the deviation a the expression 



Seff = Sb + ^ J [dacrd'^a - la') dil , (7) 

where dacrd^a = (dga)'^ + (d^a)"^/ sm^ 9 and dQ = sinOdOdip. Clearly, the spectrum of the 
quadratic part in eq.© is proportional to that in eq.®- 

The described separation between the hard and the soft modes becomes ambiguous at 
i ~ mRc, where the soft part merges into the hard one. The details of this merger would be 
unimportant if the path integral over the soft spectrum were convergent. Then the whole 
calculation would be reduced (by the Appelquist-Carazzone theorem jH]) to calculating the 
path integral with the effective action in eq.(j3]), i.e. in the effective 'low-energy' theory. 
This however is the case only in (1+1) dimensions jH], while in the discussed here (2+1) 
dimensional case the path integral with the effective action (j^ diverges, and one has to 
resort to a somewhat more accurate application of the Appelquist-Carazzone theorem, which 
requires [Ti] an explicit consideration of the regularization of the effective low-energy theory. 

In order to describe the regularization of the effective low-energy theory we concentrate 
now on the notion of the parameter fi in the limiting case of e = 0. At zero e the two 
considered vacua are degenerate and there is a stable field configuration interpolating between 
them. In the three dimensional Euclidean space the interface between the vacua makes a 
two-dimensional surface with the action proportional to the area of the surface: 

Seff = ■ Area . (8) 

The presence of the wall spontaneously breaks the translational invariance. As a result 
there appears a spectrum of (Goldstone) modes, that can be described by a massless two- 
dimensional scalar propagating on the surface of the wall. The spectrum can also be readily 
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found from the effective action (jHJ by considering small deviations of the position of the 
wall from its equilibrium (flat) shape. On a flat wall the spectrum of these modes can be 
parametrized by the two-dimensional momentum {a = 1,2): Afc = k"^. A calculation of 
the partition function in the effective 'low-energy' theory at the one-loop level immediately 
runs into the problem that the integral over the Goldstone modes is divergent. In particular, 
the renormalization of /i in the effective theory is quadratically divergent: 

In the full theory (as opposed to the effective one) however, no such divergence arises, since 
the effective description fails at ~ m, where the soft spectrum merges into the hard one. 
Thus physically the ultraviolet cutoff in the integral in eq.(jni) is at the scale m, and the surface 
tension gets a one- loop quantum correction of order m^, which is the normal behavior in the 
full theory. In order to still enable a description of the low-energy modes by the effective 
action (jH} at the loop level, it is necessary to explicitly introduce a separation parameter, 
which would serve as an ultraviolet regulator for the low-energy theory, while still being 
within the applicability of the expression for the soft modes. This can be done by using the 
standard Pauli-Villars regulator fields. 

Following the Pauli-Villars procedure we introduce a set of regulator fields ipi (at least 
two are required to regularize the quadratic divergence in eq.©) with the mass parameters 
Mj. For each regulator field the spectrum of the eigenvalues is shifted up with respect to 
those of the fields in the original field model (A„) by Mf: 

Xnm=Xn + M^. (10) 

The loop with the regulator field ^pi is given the weight factor Cj subject to the condition: 

^c, = l, ^c,M2 = 0. (11) 

i i 

The regulator mass parameters Mj are assumed to be much less than the full theory mass 
scale: Mj -C m, but much larger than the inverse size of the surface of the wall: Mj ^ 
(Area)"^/^. 

Formally, the described regulator fields are introduced in the path integral Z of the 
original theory by inserting a factor of one in the form: 

Z = ^Z, (12) 
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where = Y[i{Z[ipi]Y\ The partition functions can then be spht (at least at the one- loop 
level) into the products of the soft (Zg) and hard (Zh) contributions, e.g. Z = Z^ Z^, where 
the separation between the "soft" and "hard" modes is introduced at a scale A, intermediate 
between Mj and m: Mj <^ A ^ m. Then the identity for the partition function can be 
rewritten as 

Z=j^{iZ^),Z,}. (13) 

Clearly, the first factor is the regularized partition function described by the effective action 
in eq.®. In particular the regularized (at one loop) surface tension in this effective theory 
reads as 

(2^ = ^+87' 



1 /• 



lnA;2-^Ciln(A;2 + Mf 



where = Y.i Q Mf In Mf. 

The expression in the curly braces in eq. ()13|l is the original path integral with the soft 
modes replaced by those of the regulator fields, which implies that all the modes relevant 
for calculation of the latter expression are hard in the sense that they start at least from the 
scale of regulator masses Mj. Thus at any shape of the wall the latter expression is sensitive 
only to local properties of the surface, i.e. to higher curvatures. In particular, if the wall is 
curved with a large radius R this part can produce corrections to the effective action in eq.® 
of at most the relative magnitude 0(M~^i?~^). In particular at the curvature corresponding 
to the radius Rc of the bounce such "finite wall thickness" corrections^ are not singular in 
e at e — > 0. For this reason the discussed consideration of the regularization procedure for 
a fiat wall (at e = 0) is also applicable at the intended level of accuracy for a spherical wall 
of the bounce at a small but finite e. In the latter case the effective low-energy action (jH)) 
can be replaced by the one in eq.(^ since the term with e requires no regularization in the 
effective low-energy theory. 

Thus the discussed problem of calculating the rate of the false vacuum decay is reduced 
to evaluating the contribution of the bounce configuration and of the fluctuations around 
it to the partition function determined by the effective action (jH) and regularized by the 
'soft' regulator factor {Z^)s. As already mentioned the value of the effective action (j^ on 
the saddle point configuration reproduces the exponential factor in the decay rate, while the 



^An expression for these corrections in a i/)^ model can be found in Ref.|13|. 
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pre-exponential factor can be written as 

F = fo exp I Y.{2i + 1) \j2ciln + 1) + M^RI - u^] - In + 1) - u^]] , (15) 

where the constant /o comes from the expression in curly braces in eq. ()13|) and does not 
depend on e (or, equivalently, on Rc), and the parameter u is temporarily introduced in 
order to regularize the infrared singularity at cu^ = 2 due to the translational zero modes. 
One can also notice that, formally, the upper limit in the sums over i shout be set at a 
large value L related to the previously introduced separation parameter A as L ~ ARc. 
However the overall sum in eq. ()15|) is convergent (at imax ~ Mi) due to the Pauli-Villars 
constraints (fTTjl . and the sum can be extended to infinity as shown. Furthermore, the overall 
normalization factor in the eigenvalues is not important, since the total number of modes 
is the same for the regulator fields and the original soft part of the spectrum, so that any 
common additive term cancels in the total sum. 

The sums in eq. f|15|) can be readily evaluated (up to a numerical additive constant) using 
the Euler-Maclaurin summation formula with the result reading (at u"^ close to 2) as 

~ (Mi?J-^+V3 / 1 



F = fo |2_ ;2|3/2 [-2^ Rc) , (16) 

where In M = J2i Q InMj, and the constant /o differs from /o in eq. ()15|) only by a numerical 
factor. One readily recognizes the term proportional to in the exponent as the renormal- 
ization, according to eq. lfT^ of the surface tension /x in the effective action (jD), and this term, 
together with the contribution of the hard modes (the expression in the curly brackets in 
eq. (jl3p ) replaces the lowest-order surface tension n by the one with the one- loop correction 
in the leading semiclassical exponent. 

It is the factor with a power of MRc in eq. (fTB|) which produces a nontrivial power de- 
pendence of the pre-exponential factor F on the parameter e, through the relation R^ oc e^^. 
Although one can safely set cu^ = 2 in this term, the dependence on u is shown in eq. (fT^ in 
order to illustrate the origin of the contributions to this term: the uj'^ part comes from the 
shift of the eigenvalues with respect to i{i + 1), while the extra 1/3 in the power is due to 
the discretization of the modes on a sphere. As will be discussed few lines below, a proper 
treatment of the denominator in eq. (fT?)|) singular at cu^ = 2 produces no extra dependence of 
the factor F on e. Thus after setting the power of MR^ one arrives at the formula 

(0) for the e~^/^ behavior of the pre-exponential factor in the decay rate. 



7 



The singular at cj^ = 2 behavior arises in eq.()16|) due to the space-time translational 
invariance in the probabihty of bubble nucleation. This singularity can be readily dealt with 
by either the standard consideration |H| of the integration over the translational zero modes, 
or by using the following simple regularization procedure in terms of the effective theory |9j. 
Let us temporarily slightly break the translational invariance by introducing a dependence 
of the nucleation probability on the position x of the center of the bounce described by the 
Gaussian factor exp(— ^ x^) with a small parameter ^. The total probability of the nucleation 
in a large space-time volume is then finite and is given by 



On the other hand the shift of the center of the bounce is equivalent to an amplitude of 
the partial wave with £ = 1 of the excitation a of the surface of the bounce in eq.dZ}. Thus 
the effect of the introduced infrared regularization is equivalent to adding to the action 
([7j) additional term (3^/47r) / cr^ dQ. Although, formally, this term should be added only 
for the partial wave with i = 1 where it lifts the translational modes from zero, it can be 
safely added to all modes, since for all other modes the limit ,^ — is nonsingular. As a 
result the described infrared regularization is equivalent to the shift of the 'frequency' from 

= 2 to = 2 — 3^/{27Tfi). Using this regularized expression in eq. lfTHI) one readily finds 
from comparison with eq. ()17p that the rate F (per unit space-time volume) is finite in the 
limit ^ —>■ 0, and no new dependence on e of the factor F is introduced by the infrared 
regularization (only an extra dependence on fi does arise). 

One can see from the presented here calculation that the pre-exponential factor e^^/'^ 
in the rate of false vacuum decay in (2+1) dimensions arises in fact as an analog of the 
Casimir effect on the finite surface of the bounce due to the massless spectrum of excitations 
of the surface waves (Goldstone modes) described by the effective action (jH). As previously 
mentioned, in a (1+1) dimensional case this soft spectrum completely dominates and fully 
determines the relevant path integral for calculating the rate. In the considered here (2+1) 
dimensional case the importance of the soft part of the spectrum is greatly weakened in 
comparison, but is still sufficient for deriving the dependence on the large-scale parameter 
Rc oc e^^ of the pre-exponential factor. It further looks highly unlikely that in a (3+1) 
dimensional case the soft part of the spectrum alone can be used to make any conclusions 
about the behavior of pre-exponential factor in the rate of the false vacuum decay. 

The presented here consideration, and the resulting formula in eq.© can be used essen- 




3/2 



r . 



(17) 
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tially without modification for tlie case of thermal decay of a metastable state in a first order 
phase transition in a three-dimensional system (i.e. in (3+1) dimensions) near the stability 
point, i.e. at e — > 0. However in realistic thermal systems the nucleation of the bubbles is 
usually governed by a diffusion dynamics, rather than by the conservative Hamiltonian one 
implied throughout the discussion in this paper. The pre-exponential factor in the nucleation 
rate in a diffusion dynamics case is a well known textbook material |15J going back to the 
original work of Zel'dovich[16J. 

I thank L. Okun and A. Vainshtein for helpful discussions. I also acknowledge the support 
by the Alexander von Humboldt Foundation of my visit to the University of Bonn where this 
paper was finalized. This work is supported in part by the DOE grant DE-FG02-94ER40823. 

Note added. After this paper was completed I became aware of the work [17 , where 
the behavior of the pre-exponential factor equivalent to e~^/^ in (2+1) dimensions has been 
derived in the context of membrane creation by an antisymmetric tensor field. I believe that 
the calculation in the present paper is somewhat simpler, and is more directly related to the 
generic framework of one-loop calculations in specific field-theoretical models. 
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